Structural and vibrational properties of two-dimensional Mn x O y nanolayers on Pd(100) 
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Using different experimental techniques combined with density functional based theoretical meth- 
ods we have explored the formation of interface-stabilized manganese oxide structures grown on 
Pd(100) at (sub) monolayer coverage. Amongst the multitude of phases experimentally observed we 
focus our attention on four structures which can be classified into two distinct regimes, characterized 
by different building blocks. Two oxygen-rich phases are described in terms of MnO(lll)-like O-Mn- 
O trilayers, whereas the other two have a lower oxygen content and are based on a MnO(100)-like 
monolayer structure. The excellent agreement between calculated and experimental scanning tun- 
neling microscopy images and vibrational electron energy loss spectra allows for a detailed atomic 
description of the explored models. 

PACS numbers: 68.47.Gh, 68.37.Ef, 61.05.jh, 68.43.Bc,68.49.Uv 
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I. INTRODUCTION 

Ultra-thin oxide films grown on metal surfaces with 
thickness in the nanometer range are promising systems, 
where the reduced dimensionality and the dominance of 
surface and interface effects may lead to novel physical 
and chemical properties, intrinsically different from those 
of the corresponding bulk counterparts. As such, low- 
dimensional oxide films with tailored properties are in- 
creasingly exploited in technological applications, as for 
example in micro- and nano-electronics, magnetic record- 
ing and heterogeneous catalysis. 

From a structural point of view, systems only a few 
atomic layers thick are often different from the bulk-like 
phases in terms of local coordination, long-range order- 
ing and surface symmetry. As a consequence of the in- 
teraction with the support, even terminations which are 
not favored as bulk-truncated surfaces may be stabilized, 
and depending on the growth conditions several different 
structures appear. 

For id transition metal oxides, numerous structural 
studies of interfacial phases on various substrates are 
available in the literature, such as FeO on Pt(lll) 1-2 , 
Pt(100)2 and Ru(0001)±, vanadia on Pd(lll)£&Z. and 
Rh(lll)£&±2, Ni oxide on Pd(100)iiH and Ag(100)i2ii4, 
and Mn oxide on Pt(lll)^. In the present study we 
focus on the Mn oxide phases formed at a Mn cover- 
age of up to one monolayer (ML) on Pd(100). We have 
previously characterized the high thickness regime (15- 
20 ML), showing that MnO(lOO) with bulk- like in-plane 
lattice constant is stable in a wide range of pressure 
and temperature^; however, at more reducing condi- 
tions MnO with (111) orientation can be grown. The 
MnO(lll) surface is polar, and at high thickness the films 
tend to become unstable towards faceting: as a result, 



the morphology is dominated by triangular pyramids ex- 
posing the lowest-energy (100) facets---. At the other 
extreme of oxidizing conditions, we have demonstrated 
that flat MnO (100) films can be preferentially oxidized 
to yield Mn 3 O 4 (001)^. 

Below 1 ML, a complex phase diagram of Mn ox- 
ides on Pd(100) is observed-^, composed of nine different 
interface-stabilized phases. The phase diagram as a func- 
tion of the oxygen pressure P(02) and of the O chemical 
potential fj,Oi is depicted schematically in Fig. [T] It com- 
prises: 

- two hexagonal phases (HEX-I and HEX-II), which 
are both obtained at high O2 pressures (above 5xl0 -7 
mbar) 

- a c(4x2) structure and a stripe phase described as a 
uni-axially compressed c(4x2), which are both stabilized 
at intermediate O2 pressures (around 5xl0~ 7 mbar) 

- two structures which were called chevrons (CHEV-I 
and CHEV-II) , because of their STM appearance 

- two reduced phases with complex structures, named 
waves and labyrinth 

- finally, at the most reducing conditions, a third 
hexagonal phase (HEX-III), commensurate with the 
Pd(100) substrate (x2) along one of the two < 011 > 
directions. 

Here, we aim to rationalize this extraordinary architec- 
tural flexibility of the interfacial Mn oxides on Pd(100) 
in terms of realistic structural models. On the ba- 
sis of high-resolution electron energy loss spectroscopy 
(HREELS) data we show that it is possible to distin- 
guish two different oxygen regimes: an O-rich regime, 
comprising the HEX-I and HEX-II phases and charac- 
terized by a dominant phonon loss at w 70 meV, and 
an O-poor regime encompassing all other phases, which 
display a single phonon loss at « 44 meV. Theoretical 
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FIG. 1: Schematic phase diagram of the interfacial Mn oxides, presented as a function of the oxygen pressure P(02) and of the 
oxygen chemical potential /xo— • The nominal coverage of Mn on Pd(100) is ~ 0.75 ML. 



models based on an educated guess of possible structures 
are tested and directly compared with scanning tunnel- 
ing microscopy (STM), low energy electron diffraction 
(LEED) and HREELS results. In the following, we con- 
sider in detail four phases, c(4x2), CHEV-I, HEX-I and 
HEX-II, for which the atomic scale modeling represents 
a computationally affordable task due to the relatively 
small size of the surface unit cell. Through the com- 
bined theoretical and experimental investigation, it is 
concluded that the two O-rich phases can be described in 
terms of O-Mn-0 trilayers with a MnO(lll)-like struc- 
ture, whereas the O-poor regime is based on a compressed 
MnO(100)-like model. We find that the high degree of 
flexibility of these two basic models is favored by the for- 
mation of regular arrays of vacancies, in line with previ- 
ous investigations on related systems 12,15 . Importantly, 
our description provides a natural link between the in- 
terfacial phases and the MnO(lOO) and MnO(lll) films 
of the high thickness regime. 

The paper is organized as follows: the experimen- 
tal procedures and computational details are described 
in Section III Al and III Bl respectively. In Section IIII Al 
we discuss the main experimental evidence for the exis- 
tence of two different oxygen pressure regimes; full struc- 
tural models for the c(4x2) and CHEV-I of the O-poor 
or intermediate-0 regime are presented in Section [ill Bt 
while the proposed atomic structure for the HEX-I and 
HEX-II phases of the O-rich regime is discussed in sec- 
tion [THCl Finally, section [TVl is devoted to the summary 
and conclusions. 



II. TECHNICAL ASPECTS 
A. Experimental Procedure 

The experiments have been conducted in a number of 
different ultrahigh vacuum (UHV) chambers with typical 
base pressure p = 1-2 xl0~ 10 mbar. The structure and 
morphology of the Mn oxides films were characterized by 
STM measurements performed in our home-laboratories 
in Graz, in two different custom-designed systems. One 
STM system comprises a variable-temperature STM (Ox- 



ford Instruments), a LEED optics, a cylindrical mirror 
analyzer for Auger electron spectroscopy (AES), and the 
usual facilities for crystal cleaning and physical vapor 
deposition. The other STM system is equipped with 
a variable-temperature AFM-STM Omicron apparatus 
and a LEED optics. All STM images were recorded 
in constant current mode at room temperature (300 K) 
with electrochemically etched W tips, which were cleaned 
in situ by electron bombardment. These measurements 
were complemented by spot-profile analysis low energy 
electron diffraction (SPA-LEED) experiments performed 
in an independent UHV chamber, which has been de- 
scribed in more detail elsewhere^. HREELS measure- 
ments were performed to elucidate the phonon structure 
of the oxide layers, using an ErEELS spectrometer—, 
which was operated in specular reflection geometry at 
60° incidence from the surface normal and with a pri- 
mary energy of 6.5 eV. The crystal was kept at RT dur- 
ing the measurements, and a typical resolution of 6 meV 
was ensured, as measured at the full width at half max- 
imum of the reflected primary peak. LEED provided a 
common link with the structures prepared in the various 
chambers. 

The clean Pd(100) crystal surface was prepared by the 
usual combination of 1 keV Ar + -ion sputtering and an- 
nealing at 1000 K. Further annealing for 3-5 ruin, in 
O2 atmosphere (2xl0 -7 mbar) at 570 K, followed by a 
brief flash to 1000 K in UHV, was used to remove con- 
tamination of residual atomic carbon. The preparation 
protocol for the interfacial manganese oxide layers on 
Pd(100) consists of deposition of Mn metal in UHV with 
the sample kept at room temperature (RT), followed by 
post-oxidation (PO) at higher temperature in a molec- 
ular oxygen atmosphere. This PO procedure has been 
preferred to the reactive evaporation (RE) procedure in- 
volving deposition of Mn in oxygen atmosphere, because 
it typically provides films with smoother morphology and 
containing a single dominant phase, as judged by STM 
and LEEDi£. The Mn deposition rate was monitored 
by a quartz crystal microbalance, and evaporation rates 
of 0.5 ML min -1 or less were usually employed, 1 ML 
containing 1.32xl0 15 Mn atoms/cm 2 as referred to the 
Pd(100) surface atom density. O2 pressures in the range 
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Ixl0~ 7 -5xl0~ 6 mbar— and sample temperatures in the 
range 570-770 K were used, the exact combination of 
these parameters crucially depending on the desired long- 
range-ordered phase. 

Of the four different phases considered in detail in this 
work, the HEX-I phase is obtained by post-oxidation of 
0.75 ML Mn at 670 K at the highest oxygen pressure em- 
ployed (5xl0 -6 mbar O2), while the CHEV-I structure is 
typically prepared by post-oxidizing the same amount of 
Mn at 700 K and low 2 pressure (lx 10~ 7 mbar). Tran- 
sition between these two phases, and also to the other in- 
termediate phases c(4x2) and HEX-II, is made possible 
by annealing the HEX I phase in UHV or low O2 pressure 
at 650-700 K, or by performing a second post-oxidation 
of the CHEV-I phase at increasing O2 pressure. 

B. Computational Details 

The present calculations were carried out using the 
projector-augmented-wave (PAW) 22 ' 23 based Vienna Ab 
initio Simulation Package (VASP)2i in the framework of 
standard (Kohn Sham theory) and generalized (hybrid) 
density functional theory (DFT). The choice of adopt- 
ing a method beyond standard DFT is due to the com- 
plex nature of the system to be described. Although 
the metallic character of the Pd support is correctly de- 
scribed within DFT, this is not true for the correlated 
nature of the magnetic Mn x O y overlayers. In a series of 
recent paper s 16 ' 25 ' 26 ' 27 we have provided a detailed de- 
scription of the obstacles that standard DFT faces when 
applied to manganese oxides. We have shown that the 
application of methods beyond DFT such as HSE (Heyd- 
Scuseria-Ernzerhof) is essential to correct the self-energy 
related DFT drawbacks and to provide a satisfactory de- 
scription of the basic structural, electronic, vibrational 
and magnetic properties of these systems. 

Standard and hybrid DFT differ in the construction 
of the many body exchange and correlation functional. 
In the Kohn Sham approach we have adopted the gen- 
eralized gradient approximation scheme according to 
Perdew, Burke and Ernzerhof (PBE)2£, whereas in hy- 
brid DFT (DFTh) a suitable mixing of Hartree-Fock 
(HF) exchange with the exchange expression of PBE has 
been used following the HSE03 schem o 29 ' 30 . The result- 
ing expression for the exchange (x) and correlation (c) 
energy is 

pHSE03 1 iriHF.sr.u 1 ^ pPBE.sr.u . 7-,PBE.lr,u , piPBE 
Arc = 4 ^ + tAe + &x + h c J 

in which (sr) and (lr) denote the short- and long-range 
parts of the respective exchange interactions. The pa- 
rameter /1 controls the range separation: the optimum 
distance beyond which the short-range interactions be- 
come negligible corresponds to values of (1 between 0.2 
and 0.3 A -1 . We have used (1 = 0.3 A -1 . 

The Pd(100) substrate was modeled by repeated asym- 
metric four layers thick slabs which were used as support 



for the adsorbed Mn x O y phases. Most of the calculations 
have been performed by using 2D (4x2) and (5x3) unit 
cells employing the theoretically optimized 2D lattice 
constant. The corresponding Brillouin zone integration 
was done using 6x6x1 and 6x4x1 fc-point grids, respec- 
tively, following the Monkhorst-Pack generation scheme. 
Relaxation of all atomic positions were allowed by mini- 
mizing the atomic forces but fixing the position of atoms 
in the two deepest Pd layers. Albeit the magnetic order- 
ing of the Pd(100)-supported manganese oxide overlayers 
is unknown, throughout our theoretical investigation we 
have included spin polarization which allows for a more 
accurate description of these systems. However, to re- 
duce the already high computational costs, we have lim- 
ited our inspection to the ferromagnetic (FM) alignment, 
unless otherwise stated. The vibrational properties have 
been calculated by diagonalizing the dynamical matrix of 
the Mn x O y sub-systems. STM images have been simu- 
lated by determining iso-surfaces of constant charge, thus 
following loosely the Tersoff and Hamann method 31 . 

Finally, the stability of the most relevant models was 
underpinned by evaluating the change of the generalized 
adsorption energy 7 upon different concentration of man- 
ganese and oxygen atoms, namely: 



7 — (Eslab — -Epd(lOO) — 'TlMnfl'Mn ~ n OHo) (2) 



E s iab and -Ep^ioo) are the energies of the explored model 
and of the clean Pd(100) slab, respectively, whereas riMn 
and no are the number of Mn and O atoms for a given 
model. The reference Mn energy, /iMn is set to the energy 
of fee bulk Mn atoms (7-Mn), while /xo is treated as an 
extrinsic thermodynamic variable which can be related 
approximately to the experimental conditions under con- 
sideration (no = -1-2 eV, with respect to the free O2 
dimer). We have used a very similar approach for study- 
ing the surface phase diagram of MnO(lll)2 2 -. A detailed 
description on the thermodynamic formalism has been 
presented by Renter and Scheffler— . 

The HEX-I and HEX-II phases experimentally ob- 
served, are found to be incommensurate with (or to have 
a large coincidence mesh with respect to) the Pd(100) 
substrate. Therefore, specific to this case, we have 
performed free-standing calculations (i.e. unsupported 
Mn x O y thin layers). This is a reasonable assumption be- 
cause the interaction of oxide overlayers with the metallic 
substrate is found to be relatively weak compared to the 
interactions within the oxide itself 34 . For this case the 
stability of a particular structure was evaluated by cal- 
culating the formation energy: 



E f = (E s i ab - n Mn ^Mn - no^o)/ (riMn + n Q ) (3) 
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FIG. 2: HREELS phonon spectra of the four Mn oxide sub- 
monolayer phases discussed in this paper. From the top to 
the bottom: HEX-I, HEX-II, c(4x2) and CHEV-I. 



III. RESULTS AND DISCUSSION 
A. Evidence for two sub-monolayer regimes 

As we have briefly mentioned in the introduction, the 
phase diagram of the interfacial Mn oxides on Pd(100) is 
rather complex and made up of several distinct phases, 
which have at first sight very little in common as far as 
their long-range ordering is concerned. Identifying com- 
mon architectural aspects in these self-organized oxide 
systems is a major step on the way to rationalize the 
noticeable complexity and to unravel the driving factors 
of the architectural flexibility. HREELS offers a valuable 
diagnostic tool, because it provides the phonon loss spec- 
trum, which is tightly related to the local arrangement 
and ion coordination in the smaller building blocks con- 
stituting each long-range ordered phased. In particular, 
common building blocks shared between different oxide 
phases are expected to result in similar phonon losses, 
rendering HREELS sensitive to the occurrence of com- 
mon structural features. 

The HREELS spectra for the four phases considered 
in the present work are reported in Fig. [5] in the energy 
range up to 200 meV. All spectra exhibit a clear single 
peak structure. The phonon loss is centered at around 
70 meV for the two hexagonal phases and shifts down 
to 44-45 meV for the c(4x2) and CHEV-I phases. In- 
terestingly, the phonon spectra of the CHEV-II, waves 
and HEX III phases (not shown), obtained by further 
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FIG. 3: Calculated loss peak spectra for the MnO(100) and 
MnO(lll) unsupported thin oxide layers described in the 
text. The circle sketches provide a simplified view of the 
eigenvectors connected to the dipole active modes. Filled and 
empty circles refer to Mn and O atoms, respectively. 



lowering the O chemical potential, are also characterized 
by a single peak at 43-45 meV. These findings provide 
clear indication of two distinct regimes, an O-rich regime 
comprising phases with a higher energy phonon loss (70 
meV) and an O-poor regime with a single phonon loss 
around 44 meV. Within each regime, we may surmise 
that phases with different long-range order possess simi- 
lar building blocks^. 

Based on the hexagonal symmetry of the more oxidized 
(HEX-I and HEX-II) phases and on the lattice constant 
of the HEX-I phase, which is very close to that of bulk 
MnO(lll), one might suggest that these two interfacial 
phases are related to a MnO(lll)-like structure. Con- 
versely, the c(4x2) phase and by extension the CHEV-I 
structure might be linked to MnO(100) by analogy with 
the c(4x2)-Ni 3 O 4 /Pd(100), precursor of NiO(100)ii^ 
and obtained under similar coverage conditions. To 
seek support for this assignment we have calculated the 
dipole active modes for two free standing MnO(100) and 
MnO(lll) models containing a single manganese layer 
conveniently oxidized. The calculated phonon losses are 
reported in Fig[3] along with a sketch of the correspond- 
ing phonon eigenvectors. Indeed, the resulting phonon 
frequencies, 31 meV and 70-75 meV, unambiguously dis- 
tinguish two distinct vibrational modes characteristic of 
the (100) and (111) stacking, respectively. For these spe- 
cific models PBE and HSE provide essentially the same 
description: the main MnO(100) dipole active mode at 
31 meV originates from the out-of-plane opposite dis- 
placements of Mn and O atoms laying on the same plane, 
whereas the MnO(lll) peak at about 75 meV arises from 
the antiparallel movement of adjacent Mn and O layers. 
It is therefore natural to associate the calculated (100) 
and (111) dipole modes with the experimentally observed 
O-poor (43-45 meV) and O-rich (70 meV) phonon peaks. 
Consequently, with a certain degree of confidence we con- 
clude that the square (100) and hexagonal (111) MnO 
thin oxide building units will remain relevant in the in- 
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FIG. 4: (a) Top view of the geometrical model for the Pd(100)-supported c(4x2)-Mri304 phase: red (dark gray) spheres: 
atoms; light gray spheres: Mn atoms; green (gray) small spheres: Pd atoms underneath. Dashed lines delimit the 2D unit cell, 
the full line indicates the manganese lateral displacements and circles highlight the position of the vacancies, (b) Experimental 
STM image of the c(4x2)-Mn304 phase; sample bias U = +0.5 V; tunneling current I = 0.2 nA. (c-d) High resolution simulated 
(c) and experimental (d) STM images taken at sample bias U = +1.0 V; tunneling current I =0.25 nA. The theoretical image 
has been calculated considering tunneling into empty states between and +0.5 eV. Full lines indicate the 2D unit cell and the 
zigzag Mn chain, (e) Comparison between the measured HREELS phonon value (vertical bar) and PBE (dashed line) and HSE 
(full line) predicted dipole active modes for the c(4x2)-Mn304 phase. The inset schematically depicts the atomic displacements 
giving rise to the calculated phonon peaks. Mn and O atoms are sketched as filled and empty circles, respectively. The dashed 
line circles indicate the vacancies. 

terfacial systems grown on Pd(100). Therefore they can 
be taken as good starting models in the search for realistic 
structures of the individual phases under consideration, 
namely the c(4x2) and CHEV-I MnO(100)-like struc- 
tures and the HEX-I and HEX-II MnO(lll)-like phases. 



B. MnO(100)-like structures: c(4x2) and CHEV-I 

LEED and STM experiments have shown that an in- 
terfacial phase with c(4x2) symmetry can be stabilized 
on Pd(100) under suitable conditions 18 . This super- 



structure appears to be very similar to the recently ob- 
served c(4x2)-Ni 3 04 obtained upon reactive evaporation 
of nickel on Pd(lOO) 11 ' 12 ' 37 , which has been interpreted as 
a compressed epitaxial NiO(100) monolayer with a rhom- 
bic c(4x2) array of nickel vacancies. The analogy is not 
surprising if one considers the nature of the two oxide sys- 
tems. Indeed, strong similarities exist between NiO and 
MnO in terms of structural properties (same rhombohc- 
drally distorted fee structure), magnetic ordering (same 
antiferromagnetic type II spin ordering) and electronic 
character (partially filled localized d states). Further- 
more, bulk-like MnO(100)^ and NiO(100)^ films can 
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FIG. 5: (a) STM image showing the characteristic appearance of the chevrons structures and the local coexistence of the c(4x2), 
CHEV-I and CHEV-II phases (140 AxlOO A; U = +1.0 V; I = 0.2 nA). (b,c) Real lattice, reciprocal lattice and experimental 
LEED pattern (E=110 eV) for the c(4x2) (panel (b)) and CHEV-I (panel (c)) structures. The red arrow in the right top panel 
shows how the periodicity of the CHEV-I structure can be obtained from the c(4x2) periodicity by extending one of the two 
lattice vectors to a neighboring anti-phase position. 



be both grown on a Pd(100) substrate. 

Prompted by the strong similarities between the two 
c(4x2) phases and guided by the structural analysis on 
the c(4x2)-Ni 3 O 4 /Pd(100) phaseii>±Ml we have com- 
pared the relative stabilities of a number of different mod- 
els resembling the Ni 3 04-c(4x2) geometry. The mod- 
els are constructed by placing a compressed MnO(lOO) 
monolayer on top of a four layer Pd(100) substrate and 



creating Mn vacancies with a rhombic (RH) c(4x2) sym- 
metry. The lateral compression compensates for the pos- 
itive mismatch between MnO(lOO) (a Mn o=3.09 A, PBE) 
and Pd(100) (a Pd =2.79 A, PBE). Three options have 
been considered concerning the lateral registry: (RH1) 
oxygen atoms placed on top of Pd with the manganese 
species sitting in hollow sites; (RH2) the configuration 
opposite to RH1 with manganese arranged on top of Pd; 
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FIG. 6: (a) Top view of the geometrical model for the Pd(100) supported CHEV-I Mne07 phase: red (dark gray) spheres: O 
atoms; light gray spheres: Mn atoms; green (gray) small spheres: Pd atoms underneath. Dashed lines delimit the 2D unit cell 
and circles highlight the position of the vacancies, (b) Experimental and simulated STM images of the CHEV-I phase; sample 
bias U = +1.0 V; tunneling current I =0.13 nA. The simulated STM image has been calculated considering tunneling into 
empty states between and +1.0 eV. Unit cells are indicated by full lines, (c) Comparison between the measured HREELS 
phonon value (vertical bar) and PBE (dashed line) and HSE (full line) predicted dipole active modes for the CHEV-I phase. 
The inset schematically depicts the atomic displacements giving rise to the calculated phonon peaks. Mn and O atoms are 
sketched as filled and empty circles, respectively. The dashed line circle denotes a vacancy. 



and (RH3), with both Mn and O in bridge sites. As for 
the magnetic ordering, we have explored the FM arrange- 
ment and two different antiferromagnetic configurations 
of the Mn spins. A detailed description of the calcula- 
tions will be given elsewhere^. For the present purposes 
it is essential to point out that the FM RH1 structure, 
schematically depicted in Fig. [U^a), results as the most 
favorable one, at least 200 meV more stable than all other 
considered models. The fully optimized structure con- 
sists of a corrugated MnO(100) adlayer whose average 
distance from the substrate (2.42 A) is 22% larger than 
the Pd(100) interlayer spacing. The corrugation of « 0.23 
A arises from the outward displacement of the oxygen 
atoms. In the layer plane we find a negligible displace- 



ment of the oxygen atoms from their on-top positions, 
but a small shift (0.16 A) along the [011] direction ex- 
perienced by the Mn atoms. This displacement towards 
the vacancy clearly is the result of the free space and 
charge left by the missing manganese and acts to release 
the compressive strain in the Mn-0 over layer. The re- 
arrangements of the Mn species, also highlighted by the 
full line in Fig. 0Ja), results in a zig-zag Mn-Mn chain 
along [011] which was also observed in the c(4x2)-Ni3C>4 
system. 

In Fig. HJb)-(d) we compare experimental and calcu- 
lated STM images of c(4x2)-Mn 3 04. The close resem- 
blance between the STM simulation (inset (c)) and the 
high resolution STM image (Fig. H^d)) allows us to dis- 
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cuss further atomic details of the RH1 model. Atomic 
resolution reveals that the RH1 unit cell consists of dif- 
ferent sub-units made up of bright spots with different 
contrast and black holes spanning the 2D unit cell. The 
rhombic array of dark depressions, separated by a dis- 
tance of 5.58 A, reflects the network of Mn vacancies en- 
circled in Fig. HJa), whereas the light protrusions result 
from the remaining Mn species. The strong structural 
rearrangement induced by the lattice mismatch and by 
the creation of a Mn vacancy is reflected in the STM 
spots: the single bright features can be assigned to the 
Mn species sandwiched between two vacancies, whereas 
the light rows developing along the [Oil] direction are 
interpreted as zigzagging Mn-Mn chains embedded in a 
regular oxygen array (laying w 0.23 A above). 

Further support for the proposed c(4x2)-Mn 3 04 sur- 
face structure is provided by the comparison of the calcu- 
lated phonon spectrum with the HREELS measurements 
(Fig. 0Je)). Remarkably, the single peak structure cen- 
tered at ~ 43.5 meV revealed by HREELS is very well 
reproduced by our calculations. Although standard DFT 
underestimates the phonon energy by w 7 meV, the in- 
clusion of a portion of Hartree-Fock exchange within the 
HSE method improves significantly the agreement with 
the experimentally measured value. The diagonalization 
of the dynamical matrix enables to assign this specific 
phonon mode to the anti-phase out-of-plane vibrations 
of surface O and Mn atoms, as schematically depicted 
in the inset of Fig. [3Je). This vibrational mode is very 
similar to that found for the MnO(lOO) unsupported thin 
film (Fig. |3J) thus reflecting the MnO(100)-like nature of 
the RH c(4x2)-Mn 3 4 phase. 

Despite the striking similarity between the two c(4x2) 
monolayer structures formed by Ni and Mn oxides on 
Pd(100), precursor phases of NiO(lOO) and MnO(lOO), 
an important difference can be highlighted when consid- 
ering the lattice mismatch between the two bulk oxides 
and the Pd(100) substrate. The experimental mismatch 
is in fact «7% and 14% for NiO(lOO) and MnO(lOO), 
respectively, indicating that the c(4x2) formed by Mn 
oxide is subjected to a much higher compressive stress. 
This is reflected in the experimental STM images (Fig. 
[4jb)), which show ordered regions bearing the c(4x2)- 
Mn 3 04 periodicity separated by regions of defects. The 
incorporation of extensive defect areas into the film pre- 
sumably arises from the release of the high compressive 
strain. The latter might also explain the tendency of the 
c(4x2)-Mn 3 04 structure to nucleate preferentially near 
step edges. Moreover, we have observed that a remark- 
ably well-ordered c(4x2) overlayer forms on Pd(l 1 19)22, 
where (100) terraces exposing 10 atomic rows are sepa- 
rated by a regular array of monoatomic steps. Clearly, 
the steps provide a suitable means of lateral relaxation 
and thus enable to achieve stress relief in a natural way. 
It is also worth noting that the high surface stress in 
the c(4x2) phase may be responsible for its tendency 
to transform into a stripe phase (Fig. [T]), which can be 
interpreted as a distorted c(4x2) and is obtained at com- 
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FIG. 7: (a) Real and (b) reciprocal lattice for a (lxl) 
MnO(lll)-like hexagonal structure on a Pd(100) surface. The 
lattice parameter of the overlayer is assumed to be the mea- 
sured bulk value eiMnO= 3.14 A, whereas apd=2.75 A. In (a) 
only one of the two symmetry domains is reported for clarity, 
(c) SPA-LEED two-dimensional pattern measured at E=90 
eV, for the undistorted HEX-I phase, (d) LEED pattern of 
the distorted HEX-I phase, recorded at E=96 eV. (e) Real and 
(f) reciprocal lattice associated with the distorted hexagonal 
MnO(lll)-like phase. Only one symmetry domain is shown 
in (e) for clarity, for which the lattice parameter in the [011] 
direction of the distortion is &MnO x =2.94 A. In total, four 
symmetry domains contribute to the reciprocal lattice in (f): 
two are obtained from panel (e) with aMnO at either +60° or 
-60° from &MnO x , and the other two are obtained by rotating 
aMnO by 90° relative to the substrate mesh. 



parable values of the oxygen chemical potential^. 

We now turn our attention to the CHEV-I structure 
which is adjacent to the c(4x2) in the phase diagram 
and is obtained by lowering the O2 pressure during the 
preparation. This structure can be prepared as a single 
phase by post-oxidation method, as assessed by LEED 
and STM. However, the STM images also show that it 
can coexist locally with both the c(4x2) and CHEV-II 
phases, as illustated in panel (a) of Fig. \5\ Here, it 
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appears clearly that the c(4x2) and CHEV-I structures 
transform smoothly into each other, a strong indication 
that they possess not only similar energetics but also 
common structural features. Indeed, as it can be seen by 
inspection of the real space models in Fig. [5jb) and (c), 
the CHEV-I periodicity can be naturally obtained from 
the c(4x2) by extending one lattice vector to a neighbor- 
ing anti-phase position, thus yielding a structure with 
the matrix notation reported in Fig. 02c). As discussed 
above, the periodicity of the c(4x 2) model is given by the 
Mn vacancies. Within this picture, the CHEV-I structure 
may be described in terms of the propagation of a Mn va- 
cancy as indicated by the red arrow in Fig. 03c). Accord- 
ingly, the stoichiometry becomes Mng07, i.e. closer to 
that of MnO. Compared to the c(4x2) phase, the CHEV- 
I structure exhibits a much higher degree of long range 
order, with large defect-free areas extending over typical 
distances of 200 A and interrupted by ordered domain 
boundaries, which convey the characteristic chevron- like 
appearancei£(see also Fig. 0H a ))- 

Following the above arguments and encouraged by the 
detailed description of the c(4x2) phase in terms of the 
RH1 Mna04 model, we tested the vacancy propagation 
model adopting the above mentioned MngOy structure. 
This model is constructed by removing 14% of Mn ions 
from the compressed MnO(lOO) layer as sketched in Fig. 
UJa). The geometrical optimization yields again an oxy- 
gen terminated surface with a corrugation of 0.33 A. In 
analogy with the c(4x2)-Mn304 structure we observe a 
small displacement (0.13 A) of manganese atoms towards 
the vacancy (see arrows in Fig. [BJa)). Figure [Bfb) shows 
a comparison between the experimental and simulated 
STM images. Within the limits of the experimental res- 
olution the agreement is satisfactory. The dark circular 
depressions in the image represent the position of the Mn 
vacancies and determine the 2D unit cell of this phase, 
as outlined by the full lines. The bright protrusions are 
clearly due to manganese atoms, suggesting that elec- 
tronic effects contribute predominantly to the STM to- 
pography. Two distinct bright features are observed in 
the experimental picture, one large spot in the center of 
the unit cell and weaker flecks on the edges. The sim- 
ulated STM image together with the optimized geomet- 
rical data allow for an atomically resolved identification 
of these two kinds of features. The manganese sublayer 
is itself slightly corrugated, due to a small upward shift 
(0.05 A) of the Mn atoms laying closer to the vacancy and 
aligned in the [Oil] direction. These atoms, outlined by 
filled dots in Fig. 03(a), are responsible for the two large 
bright spots in the simulated STM image and should be 
connected to the lines of bright protrusions seen in the 
experimental STM image. The remaining spots (four per 
unit cell) arise from the lower manganese atoms. In the 
experimental STM image their intensity partially mixes 
up with the topmost Mn spots thus contributing to the 
large bright feature centered in the cell and partially pro- 
viding the weaker spots along the edges. 

Finally, our DFTh calculations establish (Fig. 0H C )) 



that the proposed Mn@07 model is characterized by a 
single dipole active mode located at around 40 meV, only 
4 meV lower than the measured loss energy. This peak 
reflects the vibrational character of the parent MnO(lOO) 
not only in terms of phonon loss energy, very similar to 
the c(4x2)-Mna04 value, but also in terms of the atomic 
displacements producing this vibration. As sketched in 
the inset of Fig. [Sfc), the peak originates from the anti- 
phase out-of-plane movements of manganese and oxygen 
atoms. 



C. MnO(lll)-like structures: HEX-I and HEX-II 

In the oxygen-rich regime, the experiments show that 
the most stable phase in a broad range of pressures above 
5x 10 -7 mbar is the HEX-I phase, which may be linked to 
MnO(lll). The structure of bulk MnO (111) is made out 
of alternating O and Mn layers, having in-plane hexago- 
nal arrangement with lattice constant OMnO- 3.14 A (see 
Fig. Efa)). A perfect MnO(lll) film would give rise to 
a hexagonal reciprocal space pattern with a lattice con- 
stant 47r/3aMnO = 2.30 A -1 , which is very close to the re- 
ciprocal Pd(100) lattice parameter (27r/apd = 2.28 A" 1 ). 
The resulting LEED pattern averaged over two hexagonal 
domains rotated by 90°, which account for the different 
symmetry of substrate and overlayer, is a circular array of 
12 extra-spots superimposed on the (10) spots of the sub- 
strate, as illustrated in Fig. 02b). Such a LEED pattern 
has been indeed observed in the experiments, as demon- 
strated in Fig. 03c)- Here, additional faint spots with 
c(6x2) symmetry are also observed, which arise from the 
coincidence mesh between substrate and overlayer—. 

A similar LEED pattern has been also obtained at 
comparable coverage for oxidized Mn on Rh(100) and 
attributed to MnO (lll)^k For Mn oxide nanolayers on 
Pd(100), however, it has not been observed always repro- 
ducibly in the different experimental systems. It appears 
that this structure is favored at a slightly lower coverage 
(0.5-0.6 ML) or after repeated oxidation cycles. Thus, 
the experimental evidence suggests that this "perfect" 
(lxl) structure may have limited stability on Pd(100) 
or may be kinetically stabilized. In contrast, the LEED 
pattern shown in Fig. OO^d) was reproducibly obtained 
after oxidation at high pressure of 0.75 ML. Here, the 
characteristic elongation of the overlayer spots indicates 
a distortion of the "perfect" (lxl) structure. It consists 
of either a contraction of one lattice vector or a reduc- 
tion of the 60° angle of the 2D hexagonal cell. The first 
type of distortion explains successfully the Moire pattern 
observed in the STM images (see Fig. 02 c )) : by reduc- 
ing the lattice vector along [Oil] to i)MnO x = 2.94 A(Fig. 
O0e)) the experimentally observed Moire periodicity of 
ps 22 A along this direction is correctly reproduced^ and 
the LEED spots appear elongated or split (Fig. 00(f))- 
By inspection of Fig. [Tfa) and (e), one may infer that 
the closer matching between the overlayer and substrate 
meshes is the driving force for the distortion: the initial 
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FIG. 8: (a-b) Top view of geometrical models for unsupported (a) HEX-I MnO(lll) and distorted (b) HEX-I MnO(lll) phases: 
red (dark gray) spheres: O atoms; light gray spheres: Mn atoms. Dashed lines delimit the 2D unit cell, (c-e) Experimental (c,d) 
and simulated (e) STM images of HEX-I; Experimental images: sample bias U = +0.5 V (c), +0.6 V (d); tunneling current 
I =0.13 nA (c), 0.15 nA (d). The simulated STM image has been calculated considering tunneling into empty states between 
and +0.5 eV. (f) Comparison between the measured HREELS phonon value (vertical bar) and PBE (dashed line) and HSE 
(full line) predicted dipole active modes for the HEX-I phase. The inset schematically depicts the atomic displacements giving 
rise to the calculated phonon peaks. Mn and O atoms are sketched as filled and empty circles, respectively. 



mismatch of 14% along the [011] direction is in fact re- 
duced to (&MnO x — apd)/apd=7%. This better epitaxial 
relationship is presumably at the origin of the increased 
stability of the distorted HEX-I phase. 

As already pointed out in Section III Bl the computa- 
tional modeling of the HEX-I and HEX-II phases has to 
be done with the unsupported setup, namely considering 
only free-standing thin layers. The most natural way to 
construct such hexagonal thin Mn-0 layers is by cleaving 
MnO perpendicular to the [111] direction. In the rock- 
salt (111) stacking the smallest blocks containing 1 ML 
of Mn atoms are either a Mn-0 bilayer or an O-Mn-0 tri- 
layer. We have first studied the relative stability of these 
two configurations by minimizing the formation energy 
with respect to the 2D hexagonal lattice constant. The 



computed energies E/, plotted in Fig. [5J clearly show 
that the trilayer structure, with a formal stoichiometry 
of MnC>2, is by far more favorable in a wide portion of the 
phase diagram than the bilayer MnO structure. Only in 
the oxygen-poor regime the bilayer gains stability with 
respect to the trilayer. However, the latter conditions 
are unrealistic compared to the experimental situation 
(« 10~ 6 mbar). 

A rigid sphere sketch of the Mn02 trilayer is given in 
Fig. [5] (a). As a direct consequence of the free standing 
nature the minimized planar lattice constant, a = 2.89 A, 
is considerably smaller than the corresponding calculated 
bulk lattice constant of 3.09 A. For the same reason, the 
optimized interlayer distance, 0.95 A, is found to be 25 
% shorter than the bulk value. In a second step, in order 
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FIG. 9: Thermodynamic DFT phase stability diagram of ex- 
plored HEX-I models in equilibrium with an O particle reser- 
voir controlling the chemical potential fio- The top scale is 
converted into oxygen pressure. 



to simulate the experimentally observed distorted hexag- 
onal structure we have computed the energy required to 
distort the perfect hexagonal trilayer by shrinking the 
lattice constant b by 7%. We found that the distorted 
trilayer is only 70 meV less stable than the perfect hexag- 
onal trilayer, as shown in the phase diagram of Fig. [5] 
Therefore, although wc find that the perfect hexagonal 
structure is the most favorable, only a low energy cost is 
required to convert it into a distorted hexagonal struc- 
ture. Considering that our model neglects the effect of 
the Pd(100) substrate and that the distortion allows for 
a better lattice matching with the substrate, we conclude 
that the distorted trilayer model of the HEX-I is also sup- 
ported by our calculations. We have tested the soundness 
of the proposed Mn02 trilayer model by comparing the 
experimental STM and HREELS data with the theory. 
The results are presented in Fig. [5] (c-f). As evident 
from Fig. [5] (c) the HEX-I phase displays the Moire- 
like superstructure already discussed 1 ^. The atomically 
resolved STM shown in panel (d) provides information 
on the local hexagonal symmetry of the array of bright 
protrusions, while panel (e) gives the result of the simula- 
tion. We find a fair agreement between the two images, 
and the theoretical data allow us to assign the bright 
spots to the topmost oxygen atoms. Convincing addi- 
tional support for our model comes from the comparison 
between the experimental and calculated phonon spectra 
of Fig. M (f). DFT and DFTh locate the distinct dipole 
active mode between 72 and 75 meV, which is only a 
few meV higher than the measured value of 70.5 meV. 
The displacements corresponding to this mode, sketched 
in the inset, clearly show that this particular vibration 
is linked to the alternately stacked hexagonal layers of 
manganese and oxygen atoms along the direction per- 



pendicular to the surface and is not compatible with the 
MnO(100)-like monolayer. It is also worth noting that 
the trilayer model is consistent with the O Is core level 
photoemission data presented in Ref^, which hint at the 
presence of two O Is components separated by rougly 0.4 
eV in binding energy. In fact, the stronger photoemission 
component at 529.1 eV can be assigned to the surface O 
layer, whereas the weaker line at higher binding energy 
can be related to the O layer at the interface with the 
Pd(100) substrate, in close analogy with previous assign- 
ments for trilayer oxide system s 42 ' 43 . 

The O-Mn-0 MnO (111) trilayer model is the natu- 
ral starting point for the structural understanding of the 
HEX-II phase which is obtained from the HEX-I phase 
upon reduction of the oxygen partial pressure. We rely 
on our previous work on the reconstruction of the bulk- 
terminated MnO(lll) surface 32 , which allows us to re- 
strict significantly the structural search. The STM and 
LEED analysis provide clear evidence for a well ordered 
2D hexagonal unit cell with lattice constant a w 6.0 
A (Refi 1 ^), corresponding to a compressed (2x2) recon- 
struction of the MnO(lll) termination. Our DFT cal- 
culations revealed that all reconstructions with (2x2) 
periodicity are more stable than others with (2x1) or 
(V3x y/S) symmetry. Amongst the explored (2x2) struc- 
tures we found that an oxygen terminated octopolar re- 
construction and a so-called a phase display the lowest 
surface energies in a wide range of oxygen partial pres- 
sure. Hence, we have focused our search on related (2x2) 
structures, limiting the exploration to thin layers con- 
taining < 1ML Mn. The octopolar reconstruction con- 
tains local pyramidal structures obtained by removing 
one quarter of atoms from the sub-surface and three quar- 
ters from the surface layer. It consists of triangular based 
pyramids with either O — Mn3 — O4 or Mn — O3 — Mn4 
stacking, respectively. Instead, the a structure is ob- 
tained by removing the topmost atom from the octopolar 
reconstruction, resulting in a Mn3 — O4 (Mn-terminated) 
or O3 — Mn4 (O-terminated) stoichiometry. In addition 
to the stoichiomctrically perfect octopolar and a recon- 
structions we have also considered slightly modified mod- 
els, obtained by suitable addition or removal of one sur- 
face or sub-surface atom. In Fig. [TU] we summarize the 
phase diagram of the most stable models explored. We 
first point out that all quasi-bilayer models (with ad- 
atoms) such as HEX-II (c) and HEX-II (d) are less fa- 
vorable. Conversely, the models that maintain a trilayer 
nature compatible with the octopolar and a reconstruc- 
tions are clearly preferred, such as HEX-II (a), HEX-II 
a and HEX-II (b). 

Let us now focus on the two most stable models, HEX- 
II (a) and HEX-II a, whose schematic view is given in 
Fig. [H](a) and (b). HEX-II (a) is obtained by adding one 
manganese atom on top of the HEX-I trilayer, whereas 
HEX-II a is the stoichiometrically perfect oxygen termi- 
nated a reconstruction, built by removing one surface 
oxygen atom per (2x2) unit cell from the HEX-I struc- 
ture. At first glance HEX-II (a) appears to be the most 
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FIG. 10: Thermodynamic DFT phase stability diagram of 
explored HEX-II models in equilibrium with an O particle 
reservoir controlling the chemical potential no ■ The top scale 
is converted into oxygen pressure. 

favorable model, but a deeper inspection leads to a dif- 
ferent conclusion. We first note that the Mn coverage of 
HEX-II (a) is 1.25 ML, i.e. larger than the experimen- 
tal one, and that the optimized planar lattice constant 
for the HEX-II (a) model, 5.78 A, is smaller than the 
measured value of 6.0 A. Conversely, the O-terminated 
a model has the right Mn coverage of 1 ML, and its re- 
laxed lattice constant is 5.95 A, i.e. almost identical to 
the experimental value. Moreover, compelling evidence 
for the O-terminated a phase comes from the comparison 
between the simulated STM and HREELS data and the 
corresponding experimental results given in Fig. [11] (c) 
and (d), respectively. The STM image reveals a hexag- 
onal pattern of large bright triangular intensity maxima 
well reproduced by the simulated STM image. The latter 
display triangular protrusions, which are ascribed to the 
oxygen trimers outlined in the sketch of panel (b). In- 
stead, the simulated HEX-II (a) STM image (not shown) 
is characterized by circular bright protrusions. As for 
the HREELS spectrum, in analogy with its precursor, 
the HEX-I phase, the HEX-II a displays a single peak 
centered at 69.5 meV, characteristic of the MnO(lll)- 
like structures. The simulation reproduces well the po- 
sition of the vibrational peak and ascribes it to the con- 
certed movement of the surface oxygen trimers against 
the manganese atoms underneath. In contrast, the cal- 
culated dipole active mode for the HEX-II (a) model is 
significantly higher in energy, 83 meV, and comes from 
a different kind of displacements, namely the vertical vi- 
bration of the octopolar Mn apex against the triangular 
oxygen basis. Taken together, it appears that the a struc- 
ture is the best of the explored models in reproducing 
the experimental scenario in terms of structural proper- 
ties, vibrational spectra and STM images. Although the 



HEX-II (a) phase is energetically competitive and even 
somewhat preferred, other calculated properties such as 
the dipole active mode and STM image do not match as 
well with the experimental data. 



IV. SUMMARY AND CONCLUSIONS 

The Mn x O y /Pd(100) system possesses a complex and 
intriguing phase diagram. In previous wor k 16 ' 17 we have 
described the high coverage regime (15-20 ML), where 
well-ordered bulk-type MnO(100) films can be grown in 
a wide range of temperature and pressure, which in turn 
can be converted either into MnO(lll) or Mn3O4(001) 
structures by tuning the appropriate conditions of tem- 
perature and pressure. At (sub)monolayer coverages an 
even richer phase diagram is found, which comprises nine 
different structures (Fig. [I}. Here, by means of com- 
bined experimental and theoretical tools we have carried 
out a detailed investigation on four prototypical phases 
(c(4x2), CHEV-I, HEX-I and HEX-II), which belong to 
two distinct oxygen pressure regimes and are character- 
ized by well-defined structural and vibrational properties. 

At oxygen-poor conditions a two-dimensional c(4x2)- 
Mn304 structure is stable, which is described as a 
MnO(100) monolayer with a rhombic array of manganese 
vacancies, in close analogy to the c(4x2)-Ni3O-4/Pd(100) 
structur e 11 ' 12 ^ 7 . Upon reduction, vacancy propagation 
leads to the MngCv CHEV-I phase, also representing a 
monolayer structure. The close structural relation be- 
tween CHEV-I and c(4x2) is reflected in the common 
vibrational fingerprint, with a single HREELS loss peak 
at w 44 meV due to the anti-phase out-of-plane move- 
ments of almost co-planar Mn and O ions. 

In the oxygen-rich regime the scenario changes sub- 
stantially. The observed HEX-I and HEX-II phases pos- 
sess a MnO(lll)-derived O-Mn-0 trilayer structure and 
display the same dipole active mode at « 70 meV. The 
transition between HEX-I and HEX-II is driven by the 
formation of oxygen vacancies inducing a so-called a re- 
construction and changing the stoichiometry from Mn02 
to M^Cv. Incidentally, the occurrence of similar trilayer 
models have been also invoked to describe other transi- 
tion metal (Rh and Ir) surface oxide a 42 ' 43 i 44 . 

The correct assignment and atomic-level description 
of all structures is validated by the excellent agreement 
between experiment and theory in terms of structural 
and vibrational properties. 
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FIG. 11: (a-b) Top view of the two most stable models for unsupported HEX-II, namely (a) octopolar-like MnsOg and 
(b) Mn407-a phases: red (dark gray) spheres: O atoms; light gray spheres: Mn atoms. Dashed lines delimit the 2D unit 
cell, whereas triangles highlight the octopolar Mn-03 pyramids and the O3 trimers of the a termination. (c) Experimental and 
simulated STM images of the hexagonal a HEX-II phase; sample bias U = +0.75 V; tunneling current I =0.1 nA. The simulated 
STM images have been calculated considering tunneling into empty states between and +0.75 eV. Unit cells are indicated by 
full lines, (d) Comparison between the measured HREELS phonon value (vertical bar) and PBE (dashed line) and HSE (full 
line) predicted dipole active modes for the a HEX-II phase. The inset schematically depicts the atomic displacements giving 
rise to the calculated phonon peaks. Mn and O atoms are sketched as filled and empty circles, respectively. 



* Electronic address: cesare.franchini@univie.ac.at 
t Electronic address: francesco.allegretti@uni-graz.at 

1 M. Ritter, W. Ranke, and W. Weiss, Phys. Rev. B 57, 
7240 (1998). 

2 L. Giordano, G. Pacchioni, J. Goniakowski, N. Nilius, E. 
D. L. Rienks, and H.-J. Freund, Phys. Rev. B 76, 075416 
(2007). 

3 S. Shaikhutdinov, M. Ritter, and W. Weiss, Phys. Rev. B 
62, 7535 (2000). 

4 G. Ketteler and W. Ranke, J. Phys. Chem. B 107, 4320 
(2003). 

5 S. Surnev, L. Vitali, M. G. Ramsey, F. P. Netzer, G. 
Kresse, and J. Hafner, Phys. Rev. B 61, 13945 (2000). 



S. Surnev, G. Kresse, M. G. Ramsey, and F. P. Netzer, 
Phys. Rev. Lett. 87, 086102 (2001). 

7 S. Surnev, G. Kresse, M. Sock, M. G. Ramsey, and F. P. 
Netzer, Surf. Sci. 495, 91 (2001). 

8 J. Schoiswohl, S. Surnev, F. P. Netzer, and G. Kresse, J. 
Phys.: Condens. Matter 18, Rl (2006). 

9 J. Schoiswohl, G. Kresse, S. Surnev, M. Sock, M. G. Ram- 
sey, and F. P. Netzer, Phys. Rev. Lett. 92, 206103 (2004). 

10 J. Schoiswohl, M. Sock, S. Eck, S. Surnev, M.G. Ram- 
sey, F. P. Netzer, and G. Kresse Phys. Rev. B 69, 155403 
(2004). 

11 S. Agnoli, M. Sambi, G. Granozzi, A. Atrei, M. Caffio, and 
G. Rovida, Surf. Sci. 576, 1 (2005). 



14 



12 S. Agnoli, M. Sambi, G. Granozzi, J. Schoiswohl, S. 
Surncv, F. P. Netzer, M. Ferrero, A. M. Ferrari, and C. 
Pisani, J. Phys. Chem. B 109, 17197 (2006). 

13 M. Caffio, B. Cortigiani, G. Rovida, A. Atrei, and C. Gio- 
vanardi, J. Phys. Chem. B 108, 9919 (2004). 

14 M. Caffio, A. Atrei, B. Cortigiani, and G. Rovida, J. Phys.: 
Condens. Matter 18, 2379 (2006). 

15 C. Hagendorf, S. Sachert, B. Bochmann, K. Kostov, and 
W. Widdra, Phys. Rev. B 77, 075406 (2008). 

16 F. Allegretti, C. Franchini, V. Bayer, M. Leitner, G. Part- 
eder, B. Xu, A. Fleming, M. G. Ramsey, R. Podloucky, S. 
Surnev, and F. P. Netzer, Phys. Rev. B 75, 224120 (2007). 

17 V. Bayer, R. Podloucky, C. Franchini, F. Allegretti, B. Xu, 
G. Parteder, M. G. Ramsey, S. Surnev, and F. P. Netzer, 
Phys. Rev. B 76, 165428 (2007). 

18 F. Li, F. Allegretti, G. Parteder, C. Franchini, R. Pod- 
loucky, S. Surnev, and F. P. Netzer, submitted (2008). 

19 J. Schoiswohl, W. Zheng, S. Surnev, M. G. Ramsey, G. 
Granozzi, S. Agnoli, and F. P. Netzer, Surf. Sci. 600, 1099 
(2006). 

20 I. Kardinal, F. P. Netzer, and M. G. Ramsey, Surf. Sci. 
376, 229 (1999). 

21 Reported pressures were not corrected for the ion gauge 
sensitivity. 

22 P. E. Blochl, Phys. Rev. B 50, 17953 (1994). 

23 G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 

24 G. Kresse and J. Furthmiiller, Comp. Mat. Sc. 15, 15 
(1996). 

25 C. Franchini, V. Bayer, R. Podloucky, J. Paier, and G. 
Kresse, Phys. Rev. B 72, 045132 (2005). 

26 C. Franchini, R. Podloucky, J. Paier, M. Marsman, and G. 
Kresse, Phys. Rev. B 75, 195128 (2007). 

27 V. Bayer, C. Franchini, and R. Podloucky, Phys. Rev. B 
75, 035404 (2007). 

28 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

29 J. Heyd, E. Scuseria, and M. Ernzerhof, J. Chem. Phys. 
118, 8207 (2003). 



30 J. Paier, R. Hirschl, M. Marsman, and G. Kresse, J. Chem. 
Phys. 122, 234102 (2005). 

31 J. Tersoff and D. R. Hamann, Phys. Rev. B 31, 805 (1985). 

32 C. Franchini, V. Bayer, R. Podloucky, G. Parteder, S. 
Surnev, and F. P. Netzer, Phys. Rev. B 73, 155402 (2006). 

33 K. Reuter and M. Scheffler, Phys. Rev. B 65, 035406 
(2001). 

34 G. Kresse, S. Surnev, M. G. Ramsey, and F. P. Netzer, 
Surf. Sci. 492, 329 (2001). 

35 S. Surncv, M. Sock, G. Kresse, J. N. Andersen, M. G. 
Ramsey, and F. P. Netzer, J. Phys. Chem. B 107, 4777 
(2003). 

36 We found further evidence for the existence of two differ- 
ent regimes in X-ray photoemission and absorption exper- 
iments with synchrotron radiation. 

37 A. M. Ferrari, M. Ferrero, and C. Pisani, J. Phys. Chem 
B 110, 7918 (2006). 

38 C. Franchini, F. Allegretti, F. Li, S. Surnev, F. P. Netzer, 
and R. Podloucky, (2008), in preparation. 

39 F. Li, F. Allegretti, S. Surnev, and F. P. Netzer, to be 
published. 

40 A c(6x2) coincidence net is only obtained when aMnO is 
expanded by w 3%. A perfect MnO(lll) structure would 
rather result in a c(8 x 2) coincidence lattice. 

41 H. Nishimura, T. Tashiro, T. Fujitani, and J. Nakamura, 
J. Vac. Sci. Technol. A 18, 1460 (2000). 

42 J. Gustafson, A. Mikkelsen, M. Borg, E. Lundgren, L. 
Kohler, G. Kresse, M. Schmid, P. Varga, J. Yuhara, X. 
Torrelles, C. Quiros, and J. N. Andersen, Phys. Rev. Lett. 
92, 126102 (2004). 

43 J. Gustafson, A. Mikkelsen, M. Borg, J. N. Andersen, E. 
Lundgren, C. Klein, W. Hofer, M. Schmid, P. Varga, L. 
Kohler, G. Kresse, N. Kasper, A. Stierle, and H. Dosch, 
Phys. Rev. B 71, 115442 (2005). 

44 Y. B. He, A. Stierle, W. X. Li, A. Farkas, N. Kasper, and 
H. Over, J. Phys. Chem 112, 11946 (2008). 



